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I. INTRODUCTION 



A Monte Carlo Study of the Rare Decay B -> XJ + £~ 
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Department of Physics, University of California at Davis, Davis, CA 95616 

(February 1, 2008) 

We study, using a monte carlo approach, the rare decay B — > X s £ + £~ including effects of the 
arbitrariness of the phase between the tp amplitudes and the perturbative amplitude, fo-quark fermi 
motion inside the B meson, and experimental smearing of lepton momenta. The fermi motion of 
b-quark inside the B meson is modeled by the ACCMM model. We found that such effects reduce 
the sensitivities of the spectra of invariant mass and forward-backward asymmetry of the lepton pair 
to new physics; especially, in the neighborhood of the ip resonances. We also estimate the sensitivity 
range of Wilson coefficients with respect to the uncertainties. 
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Rare B-meson decays are of immense interests, especially after the first discovery of the penguin decay B — > X^ by 
t^J" \ CLEO jj]. No other rare B decays have been identified yet, but they are about at the door of discovery. One of these 

■ decays that involves a lepton pair in the final state B — ► X s £ + £~ (£ = e, /i) is especially interesting. It is potentially a 
window to new physics beyond the standard model (SM) and has been studied extensively in literature. In particular, 
the spectra of invariant mass and forward-backward asymmetry of the lepton pair are shown very sensitive to different 
types of new physics Q. QCD corrections have been calculated and are under control, which enables one to predict 
reliably the decay rate and spectra || . The calculation including QCD corrections is done by an effective hamiltonian 

■ approach, which we shall describe briefly in the next section. However, a complication arises from a long distance 
0\ ' contribution of B — > X s J/%jj — ^ X s £ + £~ . To include this contribution the amplitudes of ip resonances are added to the 

perturbative amplitude in a rather phenomenological way [Q . Each of the ip amplitudes has an overall normalization 
to be determined by experiments, an arbitrary phase relative to the perturbative amplitude, Q and the ip resonance 
shape is described by a Breit-Wigner prescription. All these give rise to uncertainties in the prediction of the spectra. 
Another source of uncertainty to decay spectra comes from the fact that B meson is a bound state of a b quark 
Oh! & n d a light quark. The bound state effect can be represented by a fermi motion of the b quark, which is of order of 
Aqcd- We use the popular model of Altarelli et al. || to formulate the fermi motion. Another important smearing 
\ effect comes from the resolution in the measurement of lepton momenta. 

The objective of this report is to investigate the effects of (i) the arbitrariness of the phase between the ip resonance 
amplitudes and the perturbative amplitude, (ii) fermi motion of the b quark inside the B meson, and (iii) the smearing 
of lepton momenta on the spectra of invariant mass and forward-backward asymmetry of the lepton pair. We show 
that the sensitivity is reduced, in particular, in the regions next to tp resonances. Very often new physics can be 
parameterized by Wilson coefficients Ci(Myy). We shall investigate how much the Ci(Mw) are needed to change such 
that the spectra are distinguishable from the SM ones, under the effects of the above. 

The organization of the paper is as follows. In the next section, we describe briefly the calculation framework. In 
sec. Ill, we shall describe the smearing effects of fermi motion of b quark and lepton momentum measurements. In 
sec. IV, we estimate the sensitivities of the leptonic spectra to variation of Ci{Mw)- We conclude in sec. V. 

II. CALCULATION FRAMEWORK 

The detail description of the effective hamiltonian approach can be found in Refs. Here we present the 

highlights that are relevant to our discussions. The effective hamiltonian for B — > X s £ + £~ at a factorization scale of 
order 0(m b ) is given by HH§f 



1 Although an argument Q| of unitarity indicates that the phase should be unity, we shall, however, allow a more general phase 
in order to fully estimate this uncertainty. 
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H e $ = — '~-j=-V£Vtb 



CiQi)Oi(ji) + C 7y fa0 7l fa + C m faO% G fa + CgfaOgfa + C 10 faO w fa 



(1) 



The operators Oi can be found in Ref . ||,|| , of which the 0\ and 2 are the current-current operators and O3 — 0% 
are the QCD penguin operators. 7l and 0%g are, respectively, the magnetic penguin operators specific for b —> sj 
and b — > sg. Since Og and O10 are directly involved in the decay b — > s£ + £~, we list them here: 

2 2 
°o = TTT^^Lba hi^l , O10 = -^5 q7 ^L6 q £7^ , (2) 

where L(R) = (1 =p 7s)/2. Here we neglect the mass of the external strange quark compared to the external bottom- 
quark mass. 

The factorization in Eq.(|l|) facilitates the separation of the short-distance and long-distance parts, of which the 
short-distance parts correspond to the Wilson coefficients Cj and are calculable by perturbation while the long-distance 
parts correspond to the operator matrix elements. The physical quantities based on Eq. ([!]) should be independent 
of the factorization scale /i. The natural scale for factorization is of order m,b for the decay B — * X s £ + £~ . First, at 
the electroweak scale, say Mw, the full theory is matched onto the effective theory and the coefficients Ci(Mw) at 
the VF-mass scale are extracted in the matching process. Second, the coefficients Ci(Mw) at the VF-mass scale are 
evolved down to the bottom-mass scale using renormalization group equations. Since the operators CVs are all mixed 
under renormalization, the renormalization group equations for CVs are a set of coupled equations: 

6 fa) = Ufa M W )C{M W ) , (3) 

where Ufa Mw) is the evolution matrix and Cfa) is the vector consisting of Cifa's. The calculation of the entries of 
the evolution matrix U is nontrivial but it has been written down completely in the leading order |3| . The coefficients 
Cifa) at the scale Ofab) are given by [|| 
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QG«)=5>i"7 a< (J -1—6) , (4) 

8=1 
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C^fa = n%C ry {M w ) + I (r/M _ ^i) Csg(Mw) + C 2 (M W ) £ h iV a * , (5) 

i=l 

8 

Cf G fa =v M C SG (M w ) + C 2 (M w )Y / hv a * , (6) 

i=l 

Cffa s) = C 9 farj(s) + Y(s) + Y rcs (§) , (7) 
Ciofa) = C W {M W ) , (8) 

with 77 = a s (Mw)/oisfa> s = and q 2 is the invariant mass squared of the lepton pair. The a,'s, kji's, h^s, 

and hi's can be found in Ref. ||. The functions Ci(Mw) (i = 1 — 10), Cgfa, w(s), and Y(s) can be found in Refs. 
. Here we pay more attention to the term y r os(s), which is the contribution from the cc resonances: 

Y (s)--K^ Y Tfa^tH-)M^ 



a 



i=V»(lS),...,V(6S) 



I - m^, -t- iivi^l ^. 



where = 3Ci + C 2 + 3C3 + C4 + 3C5 + Cg and the parameter K is set at \K\ = 2.3 ||J^] with a phase k, and we 
vary the phase k to allow for the uncertainty in adding this long-distance contribution to the perturbative amplitude. 
In subsequent discussions, we only include the ip(lS) and ip(2S) resonances in our calculation for simplicity. In Eq. 
@, the resonance shape is described by a scale- independent Breit- Wagner prescription. We have verified that if the 
width r^, in Eq. (JsJ) is replaced by a g 2 -dependent width T^(q 2 ) = T l j J (q 2 /M 2 ) there are no visible changes to our 
results, because the width is very narrow. 

With the hamiltonian we can write down the decay amplitude for b — > s£~£ + and the spin-averaged square of the 
amplitude is given by 



\A\ 2 = ^L. \V*V tb f 
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Since the decay width 



where the momenta of the particles arc labelled by the particle symbols and q = 
depends on the fifth power of nib, a small uncertainty in nib will create a large uncertainty in the decay rate, therefore, 
the decay rate is often normalized by the experimental semi-leptonic width: 
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dT(b -► s£ + £~) 



T(b -► X c £v) 



ds 
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f(m c /m b ) 



(\Cf\ 2 + |C 10 | 2 ) (1 + 2s) + A\Cf\ 2 —^- + 12ft (cfCl 



icff' 



(11) 



where f(z) = 1 - 8z 2 + 8z e 



24z 4 lnz. The lepton forward-backward asymmetry A is defined by the angle 



between the £ and the 6-quark in the center-of-mass frame of the lepton pair: 



dA _ ^dcos9§-J°_ 1 dco S 9§ 
ds 



3[sK(CfC* w ) + 2K(CfC* w )] 



J dcos9^j + J , dcosl 



(\Cf\ 



Id 



(1 + 2s) + A\Cf\ 22 ^ + X2K (CfCf*) 



(12) 



One comment on these spectra under Lorentz transformation is given as follows. Since in the following we shall 
include the fermi motion of the b quark, we have to boost the spectra back to the B meson rest frame. In addition, 
the B meson is not at rest in the laboratory frame, the spectra may as well be boosted to the laboratory frame. Since 
the quantity s in Eqs. ( p"l| ) and (|l^) is Lorentz- invariant, the spectrum in Eq. ( |ll"| ) is almost stable under Lorentz 
boost. The slight change is due to the fact that different fermi-motion momentum p will give different m^,, which in 
turns affects the invariant mass spectrum. On the other hand, the forward-backward asymmetry is defined solely in 
the rest frame of the lepton pair because the angle between the lepton £~ and the b quark is not Lorentz-invariant. 



III. SMEARING 



The free quark model that treats the decay of a B meson as a free b quark is an idealistic model in the sense that it 
is only correct if mj, is infinitely heavy. It has been proved by using heavy quark effective theory that the correction 
is of order 1/m^; in particular the lepton spectrum starts with n — 2 poj . This behavior can also be understood in 
terms of fermi motion (FM) of the b quark inside the B meson, which causes a small offshellness of the b quark. The 
FM model, often called ACCMM model ||, is characterized by two parameters pf and the spectator quark mass m sp . 
The b quark is assumed to have a small momentum p, which follows a gaussian distribution with the parameter pp\ 

m = ~m exp &) ' (13) 

and a normalization dp p 2 (f>(p) = 1. Energy- momentum conservation requires the b quark mass to be dependent 
on p: 

m l(p) =m% + m 2 p - 2m B ^p 2 + m 2 p . (14) 

As a consequence of this continuous range of m&, spectra will be smeared. However, the invariant mass spectrum of 
the lepton pair will be affected minimally because the invariant mass is a Lorentz-invariant quantity and the lepton 
spectrum only receives corrections of order 1/m 2 . The invariant mass spectrum with the effect of FM smearing is 
given by 

dT(B -> X s £+r) = / dpp 2 4>(p) — , (15) 
Jo " s m b 

and a similar expression is valid for the lepton forward-backward asymmetry. The parameter p max is chosen such that 
the minimum for mb(p) in Eq. ( |I4| ) goes to M^. In our numerical calculation, we use a set of values for pp and m sp , 
which is consistent with the results obtained in the spectra of B — > X s ~f, B — * X£v, and B — > X s tp ||: 
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p F = 0.54 GeV and m sp = 0.15 GeV 



(16) 



While the FM smearing is of theoretical in nature, another smearing effect comes from the measurement of lepton 
energies and momenta. This smearing effect is actually stronger than the FM smearing. Note that both the angular 
and energy measurements will be affected by detector resolution. We shall employ the following resolutions, which 
are used in the CLEO B — > X s ip measurement JllJ ) in our study: 

S -§ = + 1.9 - 0.125) % , (^) 2 = (0.0015 ft ) 2 + 0.005 2 , (17) 

where E and p t are in GeV. In our study, the 6-quark momentum p and its direction are chosen randomly and the 
above resolutions are imposed on the final state lepton momenta in a event-by-event basis. The advantage of this 
monte carlo approach is that both FM smearing and lepton momentum smearing can be combined in a event-by-event 
basis. 

The smearing effects of FM and lepton momentum measurements are demonstrated in Figs. |l|(a) and (b). It is 
clear that the regions around the resonance peaks are smeared quite significantly. In Figs. |2|(a) and (b), we show 
the effect of varying the phase k between the tp amplitudes and the perturbative amplitude. We show the results for 
K = ±l,±i, (1 ± i)/y/2. For simplicity we treat the phases for i/)(1S) and ijj(2S) to be the same. One can see that 
k = ±1 allows for maximal interference with the perturbative amplitude. Anywhere in between is possible. We then 
treat the region roughly bounded by n — ±1 curves as the uncertainty in prediction in our following discussions. 
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FIG. 1. The distributions (a) (l/r b ^x c i„)dT(B -» X s £ + £~)/d<f and (b) dA/dq 2 showing the effects of fermi motion (FM) 
of the b quark inside B meson and of the leptonic smearing. 
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The distributions (a) {l/Tt->x c iv)dT(B — > X s ^ + £ )/dg 2 and (b) dA/dq 2 showing the effects of varying the phase 
the ip and perturbative amplitudes. 



IV. SENSITIVITY TO NEW PHYSICS 

We use the two-Higgs-doublet-model II (2HDMII) and a model-independent method to illustrate the sensitivities. 
We start with the 2HDMII, which is a popular extension of the SM and provides a minimal Higgs sector for super- 
symmetry. The extra contributions to the Wilson coefficients C, depend on the charged Higgs boson mass and tan/3. 
For clarity we list the coefficients Cj-io(M\y) here |p|,p]Jl^] : 



where 



C 7 (M W ) 
C S (M W ) 
C 9 (M W ) 
C W (M W ) 



A{x t ) A(x H ) 



2 6 tan 2 (3 
D(x t ) D{x H ) 
2 6 tan 2 (3 



B(x H ) 
E{x H ) 



Y 



sin 2 9 W 
Y = C(x t ) - F(x t ) 

Z = C{x t ) H — 



Xi 



8 tan 2 /? 
x t 



8 tan 2 13 



h{x H ) 



U{xh) 



A{x) = x 
B{x) = x 
D(x) = x 
E(x) = x 
C(x) = | 



8a; 2 



5a; -7 (3a; 2 - 2x) logx 



_ 12(a;-l) 3 

5x — 3 
12(x- l) 2 ~ 
a; 2 - 5.x - 2 



2(a;-l) 4 
(3a; — 2) log a; 



_ 4(a;-l) 3 
x — 3 



6(x- l) 3 
3a; log x 

2(x- l) 4 

log a; 



4(a;-l) 2 2(a;-l) 3 _ 
a; — 6 (3a; + 2) log a; 



at- 1 



(a;-!) 2 



1 



72 tan 2 /3 



h{xu) 



(18) 
(19) 
(20) 
(21) 
(22) 

(23) 

(24) 
(25) 
(26) 
(27) 
(28) 
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F(x) 
G(x) 



x x log X 

l-x + (x-l) 2 _ 

19x 3 + 25x 2 x 2 (5x 2 - 2x - 6) log x 



\ogx 



36(a;-l) 3 
xlogx 



18(x- l) 4 



fe(x) 



l-x (l-x) 2 

38a; - 79x 2 + 47x 3 (Ax - 6x 2 + 3x 4 ) log x 



6(1 -xf (l-a-) 4 
x t = m 2 /M^ ,x H = m 2 /m 2 H± . 



(29) 
(30) 
(31) 

(32) 
(33) 



We are now ready to investigate the effects of extra charged Higgs contributions to Ci(Mw) and in turns to the 
spectra of leptonic invariant mass and forward-backward asymmetry. Since 2HDMII always pushes C-?(M\y) more 
negative, the absolute value of increases and so does the rate of b — ► sj. Using the experimental rate from CLEO: 
1 x 10~ 4 < BR(b — > 57) < 4.2 x 1CT 4 at 95%CL level 0, we limit the range of charged Higgs mass to be m H ± > 400 
GeV for all tan/3 > 1. In Fig. ||, we show the invariant mass and forward-backward asymmetry for the 2HDMII 
with m,H = 400 — 800 GeV in an increment of 100 GeV and tan /3 = 1, 10, 20, 30, and 40. Here we do not include the 
effects of fermi motion nor the leptonic smearing. It is clear that the results implied by various charged Higgs mass 
m H± ?i 400 GeV cannot be easily distinguished from the SM. 




q 2 (GeV 2 ) q 2 (GeV 2 ) 

FIG. 3. (a) Invariant mass and (b) forward-backward asymmetry for the two-higgs-doublet model II with m H ± — 400 — 800 
GeV with an increment of 100 GeV and tan = 1 — 40. The two curves of the SM with « = ±1 are also shown. 

Next, we are going to use a model-independent method by varying Ci(Mw),i = 7,9,10, hoping that it can cover 
a wide variety of models of new physics. We look at each of Ct(Mw), Cg(Mw), and C\o(Mw) while keeping the 
others at the SM value. We shall estimate the range beyond which the resulting spectra are distinguishable from 
the SM ones. We do not look at C%(Mw) because C| ff (/i) do not enter Eqs. ( |TT| ) and ( [l2| ) directly and therefore 
the limit on the range of Cg(Mw) is rather loose. We found that changing the values of (77 ) g l io(Mv^) only changes 
the normalization of the continuum part of the invariant mass spectrum, which is not easy to identify given the 
experimental uncertainties. Therefore, we concentrate on the forward-backward asymmetry. We confirm that the 
forward-backward asymmetry is more sensitive than the invariant mass spectrum to changes in C7 q(Mw). However, 
for Cio(Mw) invariant mass spectrum appears to change more than the forward-backward asymmetry, but still only 
the normalization of the continuum changes. We shall discuss it in a moment. 

First we look at Cj(Mw)- We found that the forward-backward asymmetry is rather sensitive to Cr(Mw) at the 
small q 2 region. In Fig. ^, we show curves for Cl M (M w ) re -0.2, C 7 (M W ) = 0, and C 7 (M W ) = -0.4 with k = ±1. 
The region bounded by the SM curves of n = ±1 shows more or less the uncertainty in prediction. We define a 
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Cr(Mw) is distinguishable from the SM prediction when it has a significant region not overlapping with the SM 
region. As seen in Fig. ||, both C 7 (M\y) = 0, —0.4 have a region outside the SM region. For 

C 7 {M W ) < -0.4 , or C 7 {M W ) > (34) 

the forward-backward asymmetry is further distinguishable from the SM one. However, one has to be careful about 
the range of C 7 {Mw) shown in Eq. We can apply the SM evolution to evaluate the corresponding range in 

Cf(n = nth) and we obtain Cf(m b ) < -0.45 or Cf (m b ) > -0.18, respectively. The first range Cf(m b ) < -0.45 is 
already inconsistent with the experimental rate of b — > sj (the allowed range of \C^ s (mb)\ ~ 0.2 — 0.4.) The second 
range, on the other hand, has some overlaps with the experimentally allowed range. 
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FIG. 4. The forward-backward asymmetry predicted for C 7 {M W ), C 7 (M W ) = 0, and C 7 (M W ) = -0.4 with k = ±1. The 
upper curve of each set is for « = — 1 while the lower has k = 1. 

We now look at Cg(Myy). The SM value for Cg M (Mw) ss 1.6. Using the forward-backward asymmetry we found 
that 

Cg(M w ) < , or C 9 (M W ) > 4 , (35) 

is needed in order that the resulting spectra is sufficiently different from the SM curves, as shown in Fig. || In 
the SM, C W (M W ) fa —4.5. Since the forward-backward asymmetry dA/dq 2 is roughly proportional to C W (M W ), 
as indicated by the numerator of Eq. (jl^), therefore the asymmetry will not change significantly unless Cio(Mw) 
changes sign. We found that we need a rather large change in Cxq{Mw) in order for the forward-backward asymmetry 
to be distinguishable from the SM curves. We found, as shown in Fig. |^, 

C w (M w )<-8, or C 10 (M W )>-1 (36) 

is needed. However, this difference is only marginal and only at the large q 2 region, where the event rate is relatively 
low. 

Overall, we have found that we need a rather large change in Cg ) io(Mv^) in order to make the forward-backward 
asymmetry distinguishable from the SM prediction. Although C 7 (Mw) does not need to change a lot for the effect 
to be seen, the sensitivity range is, however, severely limited by the experimental rate of b — > sj. 
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FIG. 5. The forward-backward asymmetry predicted for Cg M (M w ), Cg(M w ) = 0, and Cg(Mw) = 4 with k = ±1. The 
upper curve of each set is for « = — 1 while the lower has k = 1. 
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FIG. 6. The forward-backward asymmetry predicted for Cf^lMw), Cio(Mw) = —1, and Cio(Mw) = —8 with k = ±1. 
The upper curve of each set is for k = — 1 while the lower has k = 1. 
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V. CONCLUSIONS 



In this report, we have studied the sensitivities of invariant mass and forward-backward asymmetry of the lepton 
pair in the decay B — ► X s £ + £~ to changes in the Wilson coefficients C7 9 io(Mjy)j under both the theoretical uncer- 
tainties, including the effect of 6-quark fermi motion inside B meson and the unknown phase between the perturbative 
amplitude and the long-distance ip amplitudes, as well as the experimental uncertainty of the measurement of lepton 
momenta. All these uncertainties make the SM prediction become a broad "region" that only when new physics pre- 
dictions go beyond this region can one say the spectrum is sensitive to new physics. We found that the sensitivity of 
the lepton forward-backward asymmetry is rather weak to changes in Cg t io(Mw)- Only when C$,\q{Mw) change sub- 
stantially will the asymmetry be distinguishable from the SM prediction. For Cj(Mw) the lepton forward-backward 
asymmetry is more sensitive but, however, a large part of the sensitivity range is already ruled out by the b — ■> 
rate. 
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